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Abstract (cont'd) 


(dynamic pressures of 12.7 and 17.8 cni. of H 2 ° ) > v^hich gives a Re^ range of 
2.93 x 10^ to 3.38 x 10 6 . For the Law of the Wall, Defect Law, and the turbulence 
quantities, very good agreement was found between the present results and 
those from well-established studies for a smooth, solid surface. The sand- 
paper-roughened solid wall and solid rough bonded screening wall tests showed 
a 20 ^ 30% increment in local skin friction and a slight shift in the log 
region of the Wall Law, as well as an increase in turbulence quantities over 
the smooth wall results. These results are in accord with the existing results 
for rough solid surfaces in this range. The effect of porosity was shown 
by comparing the sintered metal, porous wall results to the sand-roughened, 
solid wall results. Although there is a difference in roughness patterns 
for these two cases, the average k is in the same range of 5 * 7. To check 
the effect of porosity directly without any confusion from different surface 
roughness patterns, one can compare the results between the "smooth" per- 
forated titanium wall and the smooth, solid wall, or between the porous 
and solid bonded screening walls. The effect of porosity showed a 30 ^ 40% 
increment in local skin friction and a marked downward shift of the logarithmic 
portion of the Wall Law, as well as an increase in turbulence quantities 
over the smooth wall results. All these results demonstrate that a rough, 
porous wall simply does not influence the boundary layer in the same way as 
a rough, solid wall. Therefore, turbulent transport models for boundary 
layers over porous surfaces either with or without injection or suction must 
include both surface roughness and porosity effects. 
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TURBULENT BOUNDARY LAYER OVER SOLID 
AND POROUS SURFACES WITH SHALL ROUGHNESS 


By Fred Y, Kong, Joseph A. Schetz and Fayette Collier 
Virginia Polytechnic Institute and State University 


INTRODUCTION 


Over the years, considerable effort has been devoted to the study 
of turbulent boundary layer flow over porous or permeable surfaces with 
injection or suction. This work is important from the point of view of 
several practical applications such as: drag reduction on aerodynamic sur- 
faces by maintaining laminar flow at high Reynolds numbers by suction, cooling 
of turbine blades by injection, boundary layer separation control by suction, 
aircraft thermal protection, engine noise abatement, etc. Some situations 
involve induced injection or suction by condensation or evaporation, such 
as: atmospheric flows over forests, lakes, and the ocean, chemical processing, 
drying operations. Therefore, a detailed understanding of the nature of 
turbulent boundary layers with injection or suction is of significant en- 
gineering importance 

The literature in this field is rich, and there have been a number of 
review articles that summarize all but the most recent work, Refs. (1) - (5). 
Some experimental and theoretical studies, Ref. (6) - (9), have been extended 
towards determining the form of the Law of the Wall, Eq, (]), suitable for 
cases with injection or suction, since that is the foundation for all turbulence 
transport models of the inner region no matter what the level of complexity. 

A controversy has arisen, and the problem can be traced directly to the 
determination of the wall shear. The accurate determination of this quantity 


n YU* 

i = A 1 °g— +B O) 

is essential, since the so-called skin friction velocity, U* = /t w /p » enters 
directly into both coordinates on the Wall Law plot. Most workers have at- 
tempted to use either the slope of the velocity profile at the wall or the 
momentum integral equation to deduce the wall shear even though those tech- 
niques had been discredited years before for flows over solid surfaces. A 
notion of the poor reliability of those procedures can be obtained in Ref. 

(10) where the results of several variations of these techniques are compared. 
The reported results obtained in these ways for the more complicated case 
of flow over porous surfaces even with no injection or suction are further 
called into question because they do not agree with each other at the same 
Re^as shown in Fig. 1. It is a historical fact that the situation in the 
inner region of turbulent boundary layers on smooth, solid surfaces only 
became codified with the advent of direct measurements of skin friction with 
the floating element balances in the early 1950' s. 
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The accuracy of the existing semi -empirical theories for flows over 
smooth, solid surfaces that have been built upon the directly measured wall 
shear data sufficient for engineering estimation and design purposes. 

However, surfaces with some degree of roughness, which is unavoidable for 
any porous wall, present difficulties in the analysis. Specifically, even 
for roughness of intermediate size (10 < k < 70), on a solid surface, it has 
been observed that each roughness pattern behaves in a characteristic way. 

The phenomenon is more complicated for the small roughness regime (k = 5 - 10) 
as in Ref. (11). In that reference, Clauser presented the variation of the 
downward shift in the logarithmic region of the Wall Law, A (U/U*) vs. k for 
different roughness patterns. Since most of the porous surfaces necessarily 
involve some roughness, this extra complication must be taken into account 
carefully. Moreover, the effects of surface porosity, which allow small in- 
and out- flows through the surface even without applying Injection or suction, 
adds further complications that are only now beginning to be understood. 

Thus, an important part of any injection or suction study must be a documentation 
of the flow over the surfaces of interest without injection or suction, so 
that the effects of surface roughness, porosity, and their combination can 
be displayed clearly. Also, the surface pattern of the porous material 
under study must be documented. Table I shows that few papers in the literature 
really concentrate on flows over porous surfaces without injection or suction 
and investigate the influences of surface roughness and porosity independently 
before considering injection or suction. This is another reason why there 
is no current theory which is considered to be reliable. It is necessary 
to understand these phenomena before analytical models can be built and 
further extension can be made to more complicated situations without confusion. 

A previous investigation at Virginia Tech (Ref. 9) for flow over a 
sintered, porous wall of apparently small roughness had indicated a strong 
increase in wall shear without injection or suction. The values were.higher 
than one might have expected for a rough, solid surface at the same k r if one 
makes some seemingly plausible assumptions about the relationship of a non- 
uniform sintered metal surface to uniform sand roughness on a solid surface. 

The measurements were made with a floating element balance. Control tests 
with a smooth, solid surface in the same model with the same instrumentation 
had yielded results in excellent agreement with established results. It was 
decided to repeat the tests using a non-porous wall with a clearly describable 
roughness in the same nominal range to clarify the situation. The points at 
issue were: 1) the influence of the porosity of the first wall and the 

possibility of small in- and out- flows and 2) the proper method of charac- 
terizing the roughness "size" of a non-uniform surface such as produced by 
sintering. For reference, an oblique- view, electron microscope photograph 
of the actual sintered, porous surface is shown in Fig. 2. The porous material 
is formed from metal powder with 70% of the particles in the range from 40 
to 70 microns, with the rest smaller than 40 microns. If one takes 40 microns 
as a representative roughness size, a 600 grit sandpaper is a good match. 
Therefore, part of the current tests were run with the model surface covered 
with 3M 600 grit silicon carbide sandpaper. An electron microscope photo- 
graph of that material is shown in Fig. 3. In this picture, we can also see 
the glue which adheres the sand particles to the paper. The test conditions 
were such that k =5-7. This is in the same range as for the sintered, 
porous surface. 
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In order to examine the Influence of the porosity directly without 
further assumption as In the foregoing cases*, a perforated titanium plate 
(see Fig. 5} was also chosen for testing. Since the pores of the perforated 
titanium plate were made by impinging electronic beams on a smooth, titanium 
plate, the actual finished plate except for the holes themselves, is rather 
smooth compared to porous surfaces obtained by other processes. Therefore, 
the results of the perforated titanium wall can be compared with those of the 
smooth, solid wall to show the effects of the flow penetration through the 
perforated plate without any interference of distributed surface roughness. 

Further, a solid, rough surface formed of bonded, fine screen called 
Dynapore and a porous, rough Bynapore wall were tested to make another check 
on the effect of porosity. The Dynapore plate (see Fig. 5) was formed by 
diffusion bonding layers of screen together. The rough, solid Dynapore was 
made by bonding a Dynapore plate with a solid steel sheet underneath and the 
rough, porous Dynapore was made simply by supporting a Dynapore plate with 
a coarsely perforated steel sheet underneath. The effects of different rough- 
ness patterns for the sandpaper-roughened wall and the Dynapore plate are 
also manifested in the results. 

Besides the Law of the Wall, the turbulence flow field of the boundary 

1 2 ? — — 

layer should also be examined (u' , v , u'v" as a minimum) to aid the 
development of higher order turbulent transport models. 

It is also worth mentioning that documenting the turbulent flow over a 
porous surface without injection or suction has practical value in its own 
right. If one is considering a laminar flow control system using suction, 
the behavior of the flow if the suction system should fail is an important 
matter. Likewise, if an attempt is also to be made to reduce skin friction 
or heat transfer by injection, one needs to know the magnitude of any initial 
increase due to porous surface roughness to be overcome for no or low injection 
when a porous surface is employed compared to the base line case of a smooth, 
solid surface. 

In summary, the major goal of the present work is to study the effects 
of small roughness and porosity on turbulent boundary layer behavior. The 
detailed documentation of turbulent boundary layer flow over porous surfaces 
without mass transfer can be useful for further extension of the semi-empirical 
theory to include the phenomenon of mass transfer. 

To achieve this goal, six cases were studied in the current experiments: 

1 ) smooth, solid wall 

2) sand-roughened, solid wall 

3) sintered metal, porous wall 

4) perforated titanium wall 

5) solid, rough bonded screening wall (Dynapore) 

6) porous, rough bonded screening wall (Dynapore) 

Measurements of the skin friction directly and the profiles of mean velocity, 
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axial and normal turbulence intensities and Reynolds stress were made at 
several axial locations and two freestream velocities, which correspond to 
dynamic pressures of 12.7 and 17.8 cm. H 2 0. (45.1 and 53.3 m/sec.) 

EXPERIMENTAL APPARATUS & INSTRUMENTATION 
2.1 WIND TUNNEL 

The testing was conducted in the Virginia Tech stability wind tunnel. 

This is a closed-circuit, continuous-flow facility with a 2m. x 2m. test 
section. The tunnel has the capability of achieving a maximum speed of 67.0 
m./sec and Reynolds numbers of the order of 3x10° per m. with a low turbulence 
factor of 1.08. The tests were run at two different freestream velocities 
corresponding to dynamic pressure?, of 12.7 and 17.8 cm. H 2 0. (45.1 and 53.3 m/sec 

The tunnel static am' total pressures were measured by using a Pitot- 
static tubg, the dynamic pressure was shown on a Barocel Electronic Manometer, 
Model 1173 , and the tunnel temperature was measured by using a digital 
thermometer, Model 1 563-83-1 15-Hk . 


2.2 MODELS 

A long axi-symmetric model with a streamlined nose and a cylindrical 
main body was used. The model is 4.04 m. long and 40.1 cm. in diameter and 
was mounted on a stand with a streamlined surface surrounding the stand 
structure (see Fig. 6). The circular cross section of the model was chosen 
to eliminate any corner effects or three dimensional flows. The large outside 
diameter of the body was chosen so that possible transverse curvature effects 
would be minimized by keeping the boundary layer thickness divided by the 
radius less than about 1:10. The actual values ranged from about 0.075 to 
0.1175. 

The model consists of four major sections: the nose, the testing area, 

the pivot and counterweight area, and the conical afterbody. 

0 33 

The 76.2 cm. long nose, shaped by the equation r = 2.6 x and blended 
into zero slope at the base, was constructed of Styrofoam and fiberglass. It 
contains pressure taps along the center line and around the circumference 
for aligning the model with the flow. The fiberglass was sanded, painted 
and waxed to a smooth, shiny finish. 

The test area was designed to handle several configurations of solid 
and porous walls. There were three configurations, which includes six cases, 
tested in the current experiment (see Table II). The first of the three 
configurations tested had a 1.22m. long, solid wall, aluminum pipe that was 
machined smooth as shown in Fig. 6. Two test cases were run with this con- 
figuration: one was the smooth, solid wall case and the other one was the 

rough, solid wall case. The rough, solid wall was formed by carefully glueing 


^Manufactured by Datametrics 
--’Manufactured by Instrulab Inc. 
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a large sheet of 3M 600 grit silicon carbide sandpaper over the smooth, solid 
wall. A thin "fairing" was formed on the base of the nose, so that there 
was no step up onto the sandpaper. This configuration had four measurement 
stations (#1, #2, #3, #4) which were located at 24.1 cm., 48.3 cm., 72.4 cm., 
and 96.5 cm. from the base of the nose (see Fig. 7). 

There was only one test case run with the second configuration, which 
is the sintered metal, porous wall case . The method of manufacturing this 
surface was to take powdered stainless steel, pour into a flat mold, and then 
heat just to the melting point so the particles fused together. The thick- 
ness of the sintered metal, porous wall is 0.3 cm. In the second configuration, 
the 1.22 m. long, solid wall testing area was replaced by a 0.61 m. long 
aluminum solid wall section, and a 0.61 cm. long, sintered, stainless-steel, 
porous wall section as shown in Fig. 8. This configuration had only one 
measurement station (#4) located at 96.5 cm. from the base of the nose (see 
Fig. 7). 

Three test cases were run with the third configuration: one was the 
perforated titanium wall + case and the others were the solid, rough Dynapore** 
wall and the porous, rough Dynapore wall cases. The thickness of the titanium 
wall is 0.066 cm., and the diameter of each hole is 0.15 mm. The spacing 
between holes is approximately 0.63 mm. The thickness of the solid and 
porous Dynapore walls are 0.145 cm. and 0.132 cm., respectively. In the 
third configuration, the sintered metal porous wall was replaced by the 
perforated titanium wall section, then the solid Dynapore wall section, then 
the porous Dynapore wall section as shown in Fig. 8. There were three measure- 
ment stations (#2, #4, #5) located at 72.4 cm., 96.5 cm., and 119.4 cm,, from 
the base of the nose (see Fig. 7). 

A typical wall static pressure distribution is presented in Fig. 9. 

It can be seen that dP/dx in the main test area is very small. The static 
pressure over the actual measurement section was determined from the Pitot- 
static probe on the boundary layer rake. The static pressure ahead of that 
was determined from wall taps. Both showed values very near the tunnel static 
pressure. We believe the apparent jump shown is not real. Some workers use 
a criteria of dc p /dx per cm. less than 0.001 for a maximum V% error in Cf. 

Our largest value was 0.0006 at sta. 5. 

2.3 INSTRUMENTATION 


A. Skin Friction Balance 


The skin friction balance used for the direct wall shear measurements 
consists of a special, highly sensitive crystal strain gauge unit,*** (see 
Fig. 10), a 25 4 cm. long moment arm which was used to increase the sensitivity. 


♦Manufactured by Pall Trinity Micro Corporation 
♦Manufactured and kindly donated by Pratt & Whitney Aircraft Group 
**80 x 700 Mesh over 80 x 80 Mesh Manufactured by Ambac, Michigan Dynamics 
Division, United Technologies 

♦♦♦Manufactured by Kistler Morse Company, Bellevue, Washington 


a floating head, and a base mount for controlling the position of the element 
head in the model. Fig. 10 shows the skin friction balance setup for the 
second configuration. Two types of element head were used for the test as 
shown in Fig. 11. For the sand-roughened wall case, the heads were also 
carefully covered by sandpaper. 

For all the test cases, the balance was calibrated in the laboratory 
and in the model just before testing. The calibration curves agreed well 
and showed excellent linearity with little or no hysteresis. In all cases, 
the head was painstakingly aligned flush with the main surface. The gap around 
the head was selected to be 0. 75% of the diameter of the circular head 
following the suggestions of Ref. 12 to minimize the errors. 

The output signa], from the strain gauge was sent through a signal 
filter, model AF-410, to cut off the frequency below 0.1 Hz, then displayed 
on a HP 7100B strip chart recorder. 

B. Pitot Rake 


The mean velocity profiles in the boundary layer were obtained with a 
rake of small Pitot tubes (see Fig. 12). The rake was approximately 3.8 cm. 
high with one static port and 24 total pressure tubes staggered at increased 
spacing with distance from the wall. The rake was mounted on the model at 
the same downstream location as the skin friction balance and slightly to 
one side to prevent damage to the balance. 

The output pressures from the Pitot rake were recorded by using a 
Scanivalve System, which consists of four components: Scanivalve model D, 

pressure transducer model 237, power supply model B24-1.2, and electro- 
mechanical counter model 177004 020. The transducer output was read on a 
HP 7100B strip chart recorder with microvolt sensitivity (+0.1% full scale). 
The pressure transducer has; the range of + 0.25 psid and was calibrated 
within + 0.15% full scale range. 

C. Hot-Wire Anemometer 


The turbulence profiles were obtained with TSI and DISA Constant 
Temperature Hot-wire Anemometer System. The Anemometer System used in this 
work consists of the following components: 

1) Constant temperature hot-wire anemometer (TSI Model 1050). This 
anemometer has a very low noise (0.007% equivalent turbulence 
intensity) and high frequency response (500 kHz); therefore, it 
is suitable for high frequency turbulence measurements. 

2) Signal linearizer (TSI Model 1052). The voltage output is 
linearized with the flow velocity by using a linearizer. The 
convenience of making turbulence measurements is increased when 
the voltage output is linear with the flow velocity. However, 
some accuracy is lost when going through any signal conditioning 
circuit. The advantage and disadvantage of using the linearizer 
are discussed in Ref. (18). 


+Manufactured by Multimetrics Industries 
^Manufactured by Scanivalve Inc. 
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3) Correlator (TSI Model 101 5C). To be able to measure the turbulence 
Intensity in two perpendicular directions and their directional 
correlation with an x-wire probe, we have to know the sum and 
difference of the output from vhe two wires. The function of the 
correlator is to perform these analog computations. 

4) Digital integrating voltmeter (TSI Model 1076). Mean quantities 
of the flow were obtained from this voltmeter* The range of its 
integrator time constant is from 0.1 to 100 seconds and the voltage 
output range is from 0.001 to 100 volts. 

5) RMS voltmeter (DISA Model 55D35). The RMS voltmeter serves to 
read the fluctuating quantities of the flow. It has a signal 
response range from 1 Hz to 400 kHz, integrator time constants from 
0.3 to 100 seconds and an accuracy of 0.5% of full scale deflection. 

A single wire, DISA 55P14 (5 pm diarn. 1.25 mm long), was used for axial 
turbulence intensity (iPTT); and an x-wire, TSI 1243 (5 urn diam., 45°, 1 mm 
s pacin g) , was used for normal turbulence intensity (v 7 " 2 *); and Reynolds stress 
(u 1 2 3 4 5 6 v r ) . An automatic traverse was used to move the hot-wire through the 
boundary layer. The hot-wire probes were all carefully calibrated before 
and during testing. 


EXPERIMENTAL PROCEDURE 

All the tests were run at two nominal speeds of 45.1 and 53.3 m/sec 
(dynamic pressure of 12.7 cm. and 17.8 cm. of H?0) to give Reynolds numbers 
of approximately 2.93 x 10° and 3.38 x 10° based on the length to the first 
station. 

Six cases were studied in this work. Each test follows the same 
procedure as described below: 

1) The skin friction balance was calibrated in the model. Data for 
skin friction were then taken. 

2) The pressure distributions were taken. 

3) The Pitot rake was then mounted slightly to one side of the floating 
head, and the mean velocity profiles were obtained. 

4) A single-wire probe, DISA 55P14, was mounted on the automatic traverse, 
and the measurements for axial turbulence intensity profiles were 
taken. 

5) The single-wire probe, was replaced by an x-Wire probe, TSI 1243, 

and tne measurements for axial turbulence intensity, normal turbulence 
intensity, and Reynolds stress profiles were obtained. 

6) For the perforated titanium wall and the porous Dynapore wall cases, 
a miniature single-wire probe, DISA 55P11 (5 pm diam. 

was mounted inside the model and about 0.159 cm. below the surface 
to check the flow penetration through the porous surfaces. 

DATA REDUCTION & ERROR ANALYSIS 

4.1 DATA REDUCTION 

A. Skin Friction 

The calibration curve of the skin friction balance is obtained by using 
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the setup as shown in Fig. 13. The thread has been straightened before it 
was used for the calibration. From the force balance, we found 


• = W 
c ” tanB - tana 


( 2 ) 


where F is the axial force acting on the floating head and W is the weight. 

w 

The calibration curves are linear and the error is less than + Z% if the 
total axial force acting on the floating head is not over 1.0 gram’which was 
true for most of our cases. 

From Eq. (3), the wall shear could be calculated. 


T A 

w A 


(3) 


where F , the total force acting on the surface of the floating head during 
the runs, is obtained by interpolation from the calibration curve and A is the 
surface area of the floating head. Then, the skin friction coefficient was 
obtained Eq. (4). 


w 


C f = 1 /2p U 2 


Finally, the skin friction velocity was calculated from Eq. (5). 


( 4 ) 


u* -AT? 


(5) 


B. Mean Velocity 

The mean velocity was obtained from the Bernoulli's equation 

P 0 = P s + l/2p U 2 (6) 

where the total pressure P and the static pressure P s were obtained from 
the Pitot rake during the same run. Any explicit effect of turbulence on 
these mean measurements was neglected. 

C. Turbulence Quantities 

The axial turbulence intensity could be obtained both from the single- 
wire probe data from Eq. (7) or from the x-wire probe data from Eq. (8). 

However, the data obtained by the single-wire probe are generally considered 

to be more accurate than those by the x-wire probe. The normal turbulence intensity 


and Reynolds stress were obtained by Sq. (9) and Eq. (10). 
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( 10 ) 


The derivations of the above equations are detailed in Ref. 13, The 
above quantities could also be normalized by the edge velocity or the skin 
friction velocity instead of the local mean velocity through some simple 
algebraic operations. 


4.2 ERROR ANALYSIS 

The following analysis consists of calculating the maximum errors in 
the measured and computed data. 

A. Vertical Distance 

The vertical distance (Y) of each individual Pitot tube of the boundary 
layer rake is measured by, an optical traveling microscope. The accuracy is 
within + 0.00054 cm. or Y + 0.35. 

The automatic traverse used for moving the hot-wire probes in the boundary 
layer can give a minimum increment of 0.1 cm. The maximum error from calibration 
is estimated to be + 1%. The maximum estimated error for the lowest position 
of the hot-wire proBe is + 0.01 cm. Therefore, the maximum combined error in 
the determination of the probe position in tfce boundary layer could be + 0.011 cm. 
which corresponds to an error of + 5% in Y/6 . 

B. Mean Velocity 

The mean velocity profile in the boundary layer is obtained by using 
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the Pitot rake. Since the pressure transducer was calibrated within + 0.15% 
of full scale deflection and the pressure was recorded on a strip chart re- 
corder with an accuracy of + 0.1% for full scale deflection* it is estimated 
that the pressure measurements are within + 0.5% of the true value. This 
includes the error due to reading the strip chart recorder because of small 
variation from the mean value due to the likely oscillations within the 
boundary layer and tunnel vibration. Therefore, the maximum error of the 
mean velocity calculated from the pressure difference is less than 1%. 

The mean velocity profiles were satisfactory except for the results 
from the Pitot tube which sits right in the body surface, Although, the 
tip of that particular Pitot tube had been flattened, as suggested by many 
researchers, the data point on the Wall Law plot was still always higher than 
it should be. Some suggested corrections for nea, wall measurements were 
applied to our results for that point; but no substantial improvement was 
obtained. Therefore, the first data point on the Wall Law plot was ignored. 

C. Skin Friction 


The most difficult part of the experiments is to obtain consistent and 
reliable skin friction results. The errors in the skin friction results are 
due to calibration, misalignment and uncertainty in reading the voltage 
output on the strip chart recorder. 

The calibration indicated that the curves were linear and re- 
peatable. Thus, the maximum error is + 2% if the force acting on the head 
is less than 1,0 gram which is true for most of our tests. 

The slight misalignment of the element head is unavoidable for a 
cylindrical surface, although the contour of the element heads were carefully 
made to match the body surface and the 3-D positioning device was also care- 
fully designed to reduce the possibilities of misalignment. Control tests 
were done by lowering and raising the element head from its neutral position 
to see how the slight misalignment will affect the results. It is estimated 
that the maximum error will not be over + 2.5% if the element head is carefully 
aligned. 

The error due to the uncertainty in reading the strip chart recorder 
was found to be less than + 1.5%. 

Combining these errors gives a resultant maximum deviation of + 6% 
in the measured value of the skin friction. The actual apparent error obtained 
for our smooth wall tests was approximately + 2% in comparison with well- 
accepted results. 

D. Turbulence Quantities 

All the turbulence quantities were measured over the center of the 
element head. The gap around the heads and the possible slight misalignment 
of the element head may cause some errors in the measurements; however, these 
errors had been minimized by limiting the gap size and carefully aligning the 
element heads to match the contour of the body surface. This error has been 
estimated to be within + 1%, 
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The major errors in the hot-wire results are due to the calibration, 
the temperature change during the test runs, and the uncertainty in reading 
the output through the voltmeters because of the oscillations of the flow 
in the boundary layer. Errors due to calibration were found to be less than 
+ 1% mainly in adjusting the span for the linearizer. Errors in reading the 
output data through the DC voltmeter and RMS voltmeter are found to be less 
than + 1% and + 0,5% respectively. However* the errors due to the temperature 
change during the test runs could be serious; therefore, the measurements 
were generally started after the tunnel was running for a short while when 
the temperature becomes stable . {± 3°C). Also, the hot-wires were 
calibrated more often if there had been a significant temperature change. 

For the single-wire measurements, the error is not high because the test runs 
could be completed in a short time and the temperature change during this 
period is small. It is estimated to be less than + 1.5% for both DC and 
RMS outputs. For x-wire measurements, the maximum error are higher and are 
estimated to be + 3% for both DC and RMS outputs. 

From the equation in section 4.1 and the errors described above, we 
can estimate the maximum errors for axial and normal turbulence intensities, 
and Reynolds stress to be roughly + 3%, + 5% and 11% respectively. These 
figures are consistent with what was found for our smooth wall results. 

RESULTS 

5.1 SMOOTH, SOLID WALL RESULTS 

The smooth, solid wall tests were run primarily to document the adequacy 
of the instrumentation and procedures. These data also serve as "baselines" 
for comparison with the main test results. 

Some mean velocity profiles obtained with the boundary layer rake are 
shown in Fig. 14a and b and 15 in terms of the conventional Wall Law plots 
and Defect Law plots (using Y/6 , since 5 can be more precisely determined). 
The good agreement with Clauser's suggestion (Eq. (11)) for the logarithmic 
portion of the turbulent boundary layer can be seen on these plots. This 
implies accurate mean velocity profile and skin friction measurements. 




+ 4.9 


OD 



2.44 lnj+ 2.5 


( 12 ) 


where <s ~ 



U e- U 

U* 


= 9.6 (1 - |) 2 


(13) 


Also, good agreement with suitably modified Clauser (Eq. (12)) and Hama 
(Eq. (13)) formulas (Ref, 14) for the inner and outer region of the turbulent 
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boundary layer #as found on the Defect Law plots. They were modified to 
correspond to <5 rather than <$. The data obtained at two different free 
stream velocities and four different axial locations showed excellent con- 
sistency. 

The local skin friction coefficients are listed in Table IV along with 
some previous results in the same apparatus. These data, except for the new 
results at Sta. M at 12.7 cm. H 2 0, agree within + 1% with the well-established 
simple law (Ref. (15)): 

C f » 0.0128 (Re 5 *) _1/6 (14) 

Typical axial turbulence Intensity plots are shown in Figs. 16 and 17. 

The data agreed very well with Klebanoff's classical results (Ref. (16)) for 
a smooth flat plate except near the outer edge of the layer. Comparisons 
not included here showed that data obtained with the x-wire were virtually 
identical to those obtained with the straight wire. 

The axial and normal turbulence intensities obtained by the x-wire are 
compared with Klebanoff's results in Fig. 18 showing excellent agreement. 

Fig. 19 shows the normal turbulence intensity profiles at three different 
axial locations. 

Fig. 20 shows that Reynolds stress profiles obtained by the x-wire are 
a little lower than Klebanoff's results. It should be noted that Kleb anoff 
did not measure wall shear independently. He obtained values of u'v' up to 
the wall by extrapolation, thus forcing a value of unity at y = 0. Fig. 21 
shows similar trends for the Reynolds stress profiles obtained at different 
locations with wall shear measured independently. 

5.2 SAND-ROUGHENED, SOLID WALL RESULTS 

These tests were run with the surface covered with 3M 600 grit silicon 
carbide sandpaper. The roughness pattern can be seen to be very uniform on 
the electron microscope photograph in Fig. 3. If one takes 1/600 in. as the 
definition of the roughness size, then the k + is about 5 for the q = 12.7 cm. 
H 2 0 test. 

Law of the Wall plots for all stations on the sand-roughened and smooth 
walls are shown in Figs. 22a, b, c and d. It can be seen that the logarithmic 
portion of the Wall Law does shift downward slightly by an amount aU 1 1.0 
at Station 1 and 2, but it does not shift at Station 3 and shifts upward by 
an amount AU + ~ 1.0 at Station 4. Previous results for uniform sand grain 
roughness (Ref. 11) indicate that the downward shift is negligible in this 
k + range. It could perhaps be concluded that; the shifts observed are within 
our error range, although they do show a consistent trend from station to 
station and test to test. This matter may well deserve further careful study. 

The skin friction coefficients listed in Table III indicate a clear 
increase of the wall shear at all stations even for the roughness corresponding 
to this low k + range. The increase is about 13% a, 21% for q = 12.7 cm. H 2 0 
test and 21% *\» 34% for q = 17.8 cm. H o 0 test. These results are also in close 
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agreement (+ 2 %) with the data for uniform sand on flat plates deduced from pipe 
flow measurements (Ref, (17)). 

Defect Law plots for sand' roughened and smooth walls are shown in Fig. 23. 

It can be seen that the velocity Defect Law is not affected by the roughness in 
this low range. 

In Fig. 24, we show some comparisons of the axial and normal turbulence in- 
tensity profiles for rough and smooth walls. The effect is largest for the normal 
turbulence intensity. It is clear that even for this low k , sand roughness has 
a measurable effect on the turbulence field throughout the boundary layer. 

The Reynolds stress profiles for rough and smooth walls are compared in the 

usual way as shown in Fig. 25(a). Since there are variations in U* between the 2 

two cases, we believe it may be more enlightening to consider the results as u'v7U , 
where U is the local mean velocity, as shown in Fig. 25(b). In both cases, the 
rough wall results are higher as expected. 

5.3 SINTERED METAL, POROUS WALL RESULTS 

The sintered metal, porous wall results, without injection or suction serve 
to display the effects of porosity and small distributed roughness that accompany 
most of the porous surfaces. The porous material used for the present study was 
made from powdered stainless steel. The particles are not spherical originally, 
and they are distorted somewhat in the sintereing process. If an average dimension 
of 40 microns is taken as the definition of the roughness size, then for the test 
at q = 12.7 cm. H 2 0, the k + is about 5. The thickness of the sintered metal, porous 
wall is 0.3 cm. 

The sintered metal, porous wall and the sand-roughned, solid wall are in the 
same nominal roughness range of k + a* 5 for the q = 12.7 cm. HoO test under the 
simple assumptions about the relationship of a non-uniform sintered metal, porous 
surface to uniform sand-roughness on a solid surface. Therefore, the comparison 
'.of these cases will indicate the effect of porosity alone. 

The Law of the Wall plot in Fig. 26 shows that the downward shift of U + for 
the sintered metal, porous wall is much greater than that for a sand-roughened, 
solid wall in the same k range. Table III shows that the skin friction coefficients 
are also increased more sharply for the sintered metal, porous wall than those for 
the sand-roughened, solid wall. It could be concluded that the porosity does con- 
tribute to the increase of skin friction and downward shift of the logarithmic portion 
of the Wall Law. This seems to be reasonable because the normal turbulence component 
will not die out on the porous surface. 

Figs. 27a and b, and 28 indicate that the axial turbulence intensity, normal 
turbulence intensity and Reynolds stress are all somewhat higher for the sintered 
metal, porous wall than those for the sand-roughened, solid wall. Of course, all 
the turbulence quantities for the two "rough" wall cases are higher than those for 
the smooth, solid wall. 

In order to compare the results between the porous rough wall and solid rough 

wall not only at the same k but also at the same k + (= — j^), a comparison was 
made between the data for the sintered metal, porous wall with a q * 12.7 cm. 
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H 2 0 at Station 4 and for the sand-roughened, solid wall with a q * 17.8 cni. 

HpO at Station 3 where the values of U* are very nearly equal, and thus k* is 
closely matched for the two cases. Fig. 29 shows that the downward shift of 
If for the sintered metal, porous wall is greater than that for the sand- 
roughened, solid wall. Figs, 30a and b and 31 indicate that all the turbulence 
quantities are higher for the sintered metal, porous wall case, 

5.4 PERFORATED TITANIUM WALL RESULTS 

Since the surface of the perforated titanium wall is, except for the 
isolated holes, rather smooth, the influence of porosity on a turbulent 
boundary layer can be seen directly by comparing the results between the 
perforated titanium wall and the smooth, solid wall. There is no need to 
make any assumption regarding roughness size and character as in the comparison 
between the sintered metal, porous wall and the sand-roughened, solid wall 
results. 


The skin friction coefficients listed in Table III show a significant 
increase for the perforated titanium wall compared to the smooth, solid wall. 

The increase is about 30%. The skin friction coefficients obtained at Station 
3 have been discarded because of experimental difficulties. It can also be 
seen from Table III that the increase of skin friction coefficients for the 
sintered metal, porous wall is about 20% over the perforated titanium wall. 

This is close to the skin friction coefficient increase for the sand-roughened, 
solid wall compared to the smooth wall. These results indicate that both small 
roughness and porosity contribute to an increase of skin friction of approximately 
20% and 30% above the smooth wall data, respectively. 


The typical Law of the Wall plot of Fig. 32 indicates that the downward 
shift of the logarithmic portion of the Wall Law is about 3.0 * 4.0 due to 
the presence of porosity on a smooth surface. Fig. 32 also shows that there 
is no appreciable further downward shift of the logarithmic region of the 
Wall Law due to the presence of the small roughness of the sintered surface. 


Fig, 33 shows that both small roughness and porosity will increase the 
axial turbulence intensity slightly. Fig. 34 shows the same trend as Fig. 33 
where the normal turbulence intensity is increased slightly by the presence 
of small roughness and porosity. The increase of axial and normal turbulence 
intensities is about 1 * 2% based on edge velocity. It was consistently found 
that the turbulence intensities decayed more rapidly across the layer for flow 
over the perforated titanium surfaces. Perhaps this is a result of the 
regular, open-spaced porosity pattern. Fig. 35 shows that Reynolds stress is 
also increased, and the increase can reach to 100% in the near wall region. 

Theoretically, the effect of porosity is due to the normal turbulence 
component penetrating the porous surface if there is no pressure gradient 
along the body surface. To check the existence of turbulence penetration 
through the surface, a very sensitive miniature single-wire probe was mounted 
beneath the perforated titanium wall to about 0.159 cm. The results in Fig. 

36 indicate that the velocity fluctuation gets higher, but the mean velocity 
remains almost zero with increasing wind tunnel speed. The slight variation 
in mean velocity could be the result of temperature change in the wind tunnel 
during testing. 


5.5 BONDED SCREENING WALL RESULTS 

The screen mesh type of roughness of the solid, rough Dynapore wall has 
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a different nature than the particle type of roughness such as the sand- 
roughened, solid wall and sintered metal, porous wall. However, the charac- 
terization of this very fine screen mesh type of roughness is just as difficult 
as the particle type of roughness. For simplicity, if we take the diameter 
of the screen wire of the Oynapore sheet, which is approximately 0.0076 cm. 
measured from Fig. J, as the representative size dimension of this rough 
surface, hen the k is approximately 9 for the q * 12.7 cm. HgO test. 

5.5.1 SOLID, BONDED SCREENING WALL RESULTS 

The "solid", Dynapore surface was made by bonding a solid sheet underneath 
a Dynapore sheet. Since the Dynapore sheet is very flimsy, it is difficult 
to cut clean holes for the floating heads. The edges of the holes and the 
nearby surface were badly damaged on the first Dynapore sheet we received. Also, 
inserting and glueing the tube around the head (see Fig. 10) tend to distort 
the surface of the cylinder. Therefore, we could not get good skin friction 
measurements with this surface except for the last station where the damage 
was minor, As an alternative, skin friction coefficients calculated from the 
Momentum Integral Method are also listed in Table III, Profiles were measured 
at additional stations. 

The effect of the screen mesh type of roughness can be shown directly 
by comparing the results between the solid, rough Dynapore wall and the 
smooth, solid wall. A Wall Law plot is shown in Fig, 37, It can be seen 
th|t the logarithmic portion of the Wall Law is shifted downward by an amount 
AU 1.0 compared to the smooth wall results. That downward shift is a 
little higher than that produced by the small sand grain roughness. 

Fig. 38 shows that the axial turbulence intensity on the solid Dynapore 
wall is higher than that of the smooth, solid wall. Also, the axial turbulence 
intensity increases more rapidly when one gets closer to the wall. 

Figs. 39 and 40 indicate that there are almost no increases in both 
the normal turbulence intensity and Reynolds stress for the solid Dynapore 
wall case. This is surprising because the surface roughness should increase 
the mixing process in the boundary layer, 

5.5.2 POROUS BONDED SCREENING WALL RESULTS 

The porous Dynapore was made by supporting a Dynapore sheet with a thin 
perforated steel plate underneath. The thickness of the perforated steel 
plate is 0.081 cm. and it has 16 small holes (d = 0.159 cm.) per square cm.. 
Although the Dynapore material is flimsy, the damage to the surface due to 
the cutting process was found to be less for this piece. It was observed that 
there is a small hill near the front edge of measuring station 4 on the surface 
of the cylindrical body. We took data both when the element head was aligned 
with the top of the hill and then aligned with the smooth part of the body 
surface. The final skin friction coefficient in Table III is the average 
value from these two measurements. 

Since the roughness on both the solid and porous Dynapore wall is 
exactly the same, the comparison of these two cases can adequately show the 
influences of porosity on the boundary layer behavior directly. 
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The skin friction data in Table III show they are about 30 * 402 higher 
for the porous, than for the solid Dynapore wall cases. That is close to the 
increase of skin friction for the sintered metal porous wall over the sand- 
roughened solid wall. Also, the increase of skin friction coefficient for 
the porous Dynapore wall over the smooth wall is about 602, which is higher 
than the corresponding increase for a sintered metal porous surface. 

The Wall Law plots of Figs. 41 and 42 show that the downward shift of 
the logarithmic portion of the Wall Law is about 4.5 for the porous Dynapore 
wall which is also a little higher than the corresponding downward shift for 
the sintered metal porous wall. Fig, 42 shows that the downward shift of 
the logarithmic region of the Wall Law for the porous Dynapore wall is much 
greater than that for the solid Dynapore wall. 

From Fig. 43, it can be seen that the axial turbulence intensity of the 
porous Dynapore wall is about 12 (based on edge velocity) higher than that 
of the solid Dynapore wall. However, the axial turbulence intensity of both 
the porous and solid Dynapore walls is higher than that of the smooth wall. 

It is also observed that the increase of the axial turbulence intensity became 
more rapid closer to the wall. 

Fig. 44 shows that the normal turbulence intensity is about 22 (based 
on edge velocity)higher for the porous Dynapore wall than that for the solid 
Dynapore wall. It also can be noted that the increase of normal turbulence 
intensity becomes more rapid closer to the wall. The Reynolds stress profiles 
in Fig. 45 indicates the same trend as the axial and normal turbulence intensity 
profiles, but the maximum increase of Reynolds stress can be as high as 1002 
(based on local mean velocity) in the near wall region. 

As we did for the perforated titanium wall, the penetration of the tur- 
bulence through the porous Dynapore was detected by a single-wire probe under 
the porous Dynapore wall. It can be seen from Fig. 46 that there is no mean 
flow but measurable turbulence passing through the porous surface. 

5.6 PRESSURE DROP THROUGH THE POROUS MATERIALS 

In order to aid in understanding the behavior observed with the three 
different porous materials used, a measure of their "porosity" was obtained 
by determining the pressure drop through each material as a function of volume 
flow rate (or average normal injection velocity). The results are given in 
Fig. No. 47 where the wide variation for the different materials can be noted. 

DISCUSSION 


The turbulent boundary layer over a smooth wall, a rough, solid wall 
and porous, rough and "smooth" walls + was studied. The goal was to investigate 
the influences of small roughness (k - 5 - 10) and porosity independently 
and their combined effects on turbulent boundary layer behavior. Here we 
have concentrated on small roughness because of its importance in practical 
engineering applications. The effect of porosity due to the existence of the 
penetration of turbulence through porous surfaces has been confirmed. 

For the smooth, solid wall case, accurate results were obtained for the 
Law of the Wall, the Defect Law, the axial and normal turbulence intensities, 
and the Reynolds stress when compared with well-established data. 
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Two solid, rough wall cases were studied in this work. One was the 
sand-roughened, solid walj. which had a sand grain type of roughness on the 
surface with an average k * 5 for the tests. The second was the solid, rough 
Dynapore wall which has a, screen mesh type of roughness on the surface. For 
both cases, the effect of small roughness is to shift the logarithmic portion 
of the Wall Law by a small to negligible amount. However, the increase in 
the skin friction values is not negligible; the Increase is about 20% * 30$ 
above the smooth wall value. It is observed that small roughness does not 
affect the Defect Law. The axial and normal turbulence intensities are in- 
creased slightly. The Increase is about 1$ based on the edge velocity, and 
the increase of Reynolds stress is somewhat higher. 

Three different porous wall cases were studied, namely a sintered metal, 
porous wall; a perforated titanium wall and a porous Dynapore wall, The effect 
of porosity can be shown when comparing the sintered metal, porous wall results 
to the sand-roughened, solid wall results. Although the character of the 
roughness for these two cases is different, the average k is in the same 
range of 5 <v 7. To see the effect of porosity directly without any interference 
of different surface roughness patterns, one can compare the results between 
the "smooth" perforated titanium wall and the smooth, solid wall, or between 
the porous and solid Dynapore walls, where the roughness patterns of these two 
surfaces are exactly the same. The comparisons reveal that the effect of 
porosity js to shift the logarithmic region of the Wall Law downward by an 
amount AU = 3^4 from the solid wall results and to increase the skin friction 
values by about 30$ * 40$. The axial and normal turbulence intensities and 
Reynolds stress are increased, and the increase becomes more rapid when one 
gets closer to the wall. 

The combined effects of small roughness and porosity can be seen by 
comparing the results between the sintered metal, porous wall and the smooth, 
solid wall or between the porous, rough Dynapore wall and the smooth, solid 
wall. The downward shift of the logarithmic region of the Wall Law and the 
increase of the skin friction value by the combined effects of small rough- 
ness and porosity is found to be roughly the sum uf their individual effects. 
However, the turbulence quantities do not show the same trend clearly, although 
all the turbulence quantities are increased by the small roughness and porosity. 

All these results demonstrate that a porous, rough wall simply does not 
influence the boundary layer in the same way as a solid, rough wall. Therefore, 
turbulent boundary layer models for cases with injection or suction must include 
surface roughness and porosity effects; however, these effects cannot simply 
be interpolated from the results for a sand-roughened, solid wall. 
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Comparison of prior experimental ( not directly measured ) results 
for skin friction on porous surfaces without injection or suction 
from Ref. (5). 
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Electron microscope photograph of 3fl 600 grit silicon carbide 
sandpaper used in the experiments. 
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Fig. 4. Electron microscope photograph of the perforated titanium 
sheet used in the experiments. 
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Axi symmetric wind tunnel model for configuration 1 
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Fig. 8. Axi symmetri c wind tunnel model for configuration 3 
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Pressure Distribution along the Body 
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Fig. 10. Skin friction balance for configuration 2 and 3. 



ELEMENT HEAD FOR CONFIGURATION 1 . 



ELEMENT HEAD FOR 


k 4,17cm. H 



4.57cm. 


11. Element heads used in the experiments. 
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Fig. 14(a). Law of the Wall plot of smooth, solid wall results, 
all stations, q=12.7 cm. Hr,0. Clauser Log Law 
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Fig. 17. Axial turbulence intensity profiles for smooth, solid wall 
tests, all stations, q=17.8 cm. H ? 0. 


WUE MERS. 



Axial and normal turbulence intensity profiles for smooth, 
solid wall tests, station 4, q=12.7 cm. HJ3. 
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Fig. 19. Normal turbulence intensity profiles for smooth, solid wall 
tests, station 2, 3, and 4, qi=12.7 cm. H„0. 



nc l n INAL PAGF IS 

0F POOR QUALITY 



<v 

UJ 



G 

cn 

UJ 


O 

i 


8 

,c 


+ 

0 


+ 

e 


o 



i n 


03 


I 

I 

O 

■hU3 


C -j, 

1 c * 
l • CE 

■ _i 


t 


UJ 




Q 

U3 

"Cv 



08 '0 


09’ 0 OP'O CZ’O 00* cP 

2 ( *n ) / . A . nr 


w . 
p 4- 
tf) O 
OJ 

P to 
P 

r— i — ■ 
r- 3 

rtl to 
5 CJ 

~o 

•r- £ 
r— P 

0- r- 

W 3 

* t3 

x: a) 

P s- 

o rrj 
O Q. 

E E 

CO o 

o 

1- 

o 

to CO 

<u X 


4- E 
o o 
s- 

Q-i"> 

to CO 
t0 r— 
OJ II 
5- cr 
P 
tO 


lO 


•=*■ !+“ 

4- 

_ c o 
- o c 

3 *r- (O 

C P JO 
<T3 <U 
UPC 
H lO 


o 

CO 


03 


39 



CvJ 

cn 


__ 

— r 

-y 

5 

6 

o 

MO 

*—• 


h— 

l— 

►—* 

g 

G 

G 

>— 

1 — 

h- 

CO 

CO 

CO 

0 

< 

+ 



40 


Reynolds stress profiles for smooth, solid wall tests 
station 2, 3, and 4, q=12.7 cm. H^O* 



SMOOTH 



Fig. 22(a). Law of the Wall plot of smooth and sand-roughened wall tests, 
station 1, q=12.7 cm. H ? 0. 




Fig. 22(b). Law of the Wall plot of smooth and sand-roughened wall tests 
station 2 , q=l 2.7 cm. H ? 0. 
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Law of the Wall plot of smooth and sand-roughened wall tests, 
station 4, q=12.7 cm. FLO. 
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Fig. 23 . Defect Law plot of smooth and sand-roughened wall tests, 
station 4, q=12.7 cm. H ? 0. 







Fig. 24. Axial and normal turbulence intensity profiles for smooth 
and sand-roughened wall tests, station'4, q=12.7 cm. H 0. 
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Fig. 25(a). Reynolds stress prufiles for smooth and sand-roughened wall 
tests, station 4, q=12.7 cm. H ? 0. 



Fig. 25(b) . Reynolds stress profiles normalized with local mean velocity 
for smooth and sand-roughened wall tests, station 4, q=12.7 
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Fig. 26. Law of the Wall plot of smooth, sand-roughened and sintered metal, 
porous wall tests, station 4, q=12.7 cm. H o 0. 
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Fig. 2a Reynolds stress profiles normalized with local mean velocity for 
smooth, sand-roughened and sintered metal, porous wall tests, 
station 4, q=12.7 cm. H ? 0. 
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roughened, sintered metal 
tests, station 4, q=12.7 
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Reynolds stress profiles normalized with local mean velocity for 
smooth, sand-roughened, sintered metal , porous, and perforated 
titanium wall tests, station 4, q=12.7 cm. HpO. 
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Fig. 36. Penetration of turbulence through the perforated titanium 
wall 
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Fig. 37. Law of the Wall plot of solid Dynapore wall tests, station 5, 
q=17.8 cm. H p 0. Compared with Clauser’s Log Law for a smooth 
wall. c 
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Fig. 38. Axial turbulence intensity profiles for smooth and solid 
Dynapore wall tests, station 4, q=17.8 cm. H ? 0. 





Normal turbulence intensity profiles for smooth and solid 
Dynapore wall tests, station 4, q=17.8 cm. H 2 0. 







Reynolds stress profiles normalized with local mean velocity 
for smooth and solid Dynapore wall tests, station 4, q=17.8 











the Wall plot of solid Dynapore and porous Dynapore wall 
station 5, q=12.7 cm. H 0 C. 





lence intensity profiles for smooth, solid Dynapore, 
Dynapore wall tests, station 4, q=12.7 cm. 
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Fig, 44. Normal turbulence intensity profiles for solid Dynapore and 
porous Dynapore wall tests, station 4, q=12.7 cm, I^O. 







Fig. 47. Pressure Drop vs. Mass Flow for the Porous Materials 
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Calculated from Momentum Integral Method! 


